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Abstract. The Boolean algebra of regular closed sets is prominent in topol- 
ogy, particularly as a dual for the Stone-Cech compactification. This algebra 
is also central for the theory of geometric computation, as a representation 
for combinatorial operations on geometric sets. However, the issue of com- 
putational approximation introduces unresolved subtleties that do not occur 
within "pure" topology. One major effort towards reconciling this mathe- 
matical theory with computational practice is our ongoing I-TANGO project. 
The acronym I-TANGO is an abbreviation for "Intersections — Topology, Ac- 
curacy and Numerics for Geometric Objects" . The long-range goals and initial 
progress of the I-TANGO team in development of computational topology are 
presented. 



1. Introduction and brief literature review 

Throughout this paper, all sets considered will be assumed to be subsets of R'^, 
with its usual topology. The Boolean algebra of regular closed sets in R'^ will be 
denoted as 7?,(R'^). Furthermore, any regular closed set considered will be assumed 
to be compact. Any surfaces and curves considered will be assumed to be compact 
2-manifolds and 1-manifolds, respectively. All neighborhoods will be assumed to 
be open subsets of R'^ . 

The theoretical role for 7?.(R'^) was introduced into geometric computing to cor- 
rect the unexpected output seen from combinatorial operations on geometric sets 
[27|. For instance, consider the two dimensional illustration shown in Figure^ The 
original operands of A and B are indicated in Figure ^a). The unexpected output 
is shown in Figure ^b), where the expected result would have been what is shown 
in Figure ^c). 

The phenomenon shown in Figure H^b) was informally described as "dangling 
edges" The formalism that was proposed to eliminate this behavior was that 
geometric combinatorial algorithms should accept only regular closed sets as input 
and then execute the Boolean operations of meet, join and complementation on 
these operands, thereby creating only regular closed sets as output [221 • The intent 
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Figure 1. Subtraction of two sets 



was to eliminate "dangling edges" and, in principle, this should have been suffi- 
cient^. However, each operand also has a geometric representation that depends 
upon the approximation methods used to compute the results. This additional 
subtlety raises issues in both theory and computation. 

For this short article, only a brief literature review will be presented. An ear- 
lier survey on topology in computer-aided geometric design j25j is recommended as 
introductory material for topologists. The texts |15l I24| discuss the integration of 
computational geometry, shape modeling and topology. The subject of computa- 
tional topology is still a nascent and emerging sub-discipline. This article focuses 
upon the authors' particular perspective in its development. The first use of the 
terminology "computational topology" appears to have been in the dissertation of 
M. Mantyla j20j. Further contemporary views can be gained from the following 
web sites 00 El- 

Additional perspective can be gained by understanding the broader context in 
which topology has already been successfully applied to computer science. We 
mention two particularly notable successes. The first is the use of non-Hausdorff 
topology by Kopperman, Meyer, Kong, Rosenfeld, Smyth and Herman in digital 
topology for computer graphics and image processing. A good overview is readily 
available ^Hl- Similarly, the work of Mislove, Reed and Roscoe on domains explores 
variants of limits for fundamental algorithmic and programming language studies, 
continuing the expressive power and the broad applicability of the language of 
topology in denotational semantics and concurrent programming. The monograph 
|26| is recommended as an introduction. 



One elegant computational representation for the combinatorial operators is to 
assign each object a symbol and then to indicate operations in a tree referencing 
those symbols. For instance, such a tree structure could be as depicted in Figure|21 

At this level of abstraction, the mathematical theory and the computational 
representation are completely consistent, and this representation became known as 
Constructive Solid Geometry (CSG). Difficulties arose in instantiating the basic 
geometric information that is represented by the operands at the leaf nodes and, 
sometimes, in computing geometric representations at the internal nodes of the 

"'^The subtraction operation between two sets, shown as A — _B in Figure^ is not specifically a 
Boolean operation. However, the use oi A — B should be understood to be conveniently shortened 
notation equivalent to the operations A/\B' , where B' represents the standard Boolean operation 
of complementation on the operand B. 



2. Theory versus computation 



COMPUTATIONAL TOPOLOGY FOR REGULAR CLOSED SETS 



3 



\ 





A 



B C 



D 



Figure 2. Tree for (A A B) V (C - D). 



tree. In CSG, the leaf nodes are restricted to a small set of specific geometric 
objects, known as primitives. A typical collection of primitives might consist of 
spheres, parallelepipeds, tori and right circular cylinders. The critical geometric 
algorithm underlying each Boolean operation is the pairwise intersection between 
the operands. 

As the boundary of each of these primitives can be represented by linear or 
quadratic polynomials, the needed intersection between each pair of primitives was 
relatively simple and numerically stable, for most cases considered, although spe- 
cific intersections could be problematical. For instance, suppose two cylinders of 
identical radius and height were created and then positioned so that the bottom 
of one cylinder was co-incident with the top of the other cylinder. This special 
case was specifically considered in most intersection algorithms and could usually 
be processed without problem. However, if one then rotated the top cylinder a 
fraction of a degree about its axis (so that the planar co-incidence remained in- 
tact) many software systems would fail to produce any output for this problem, 
sometimes even causing a catastrophic program failure. This particular problem 
became a celebrated test case and most systems developed ad hoc methods to solve 
this cylindrical intersection problem. Yet, this was just avoiding the more serious 
issue of the fragile theoretical foundations for many intersection algorithms. People 
using CSG systems became sensitive to their limitations and continued to use them 
effectively by avoiding these challenging circumstances, although the work-arounds 
were often tedious to execute. 

The imperative, largely initiated by the aerospace and automotive industries, 
to model objects using polynomials of much higher degree than quadratic created 
a movement away from CSG systems. The alternative format was to represent 
compact elements of 7^(]R^) by their boundaries, and this became known as the 
"boundary representation" approach, or "B-rep" for short. This has become the 
dominant mode today. Again, within this clean conceptual overview, the realities of 
computation pose some subtle problems. In most industrial practice, the modeling 
paradigm was further restricted so that the boundary of an object was a 2-manifold 
without boundary. However, it was difficult to create computer modeling tools that 
could globally define 2-manifolds without boundary, though there existed excellent 
tools for creating subsets of these 2-manifolds. For example, computational tools for 
creating splines were becoming prevalent. Again, in principle, if each such spline 
subset was created with its boundaries, then the subsets could be joined along 
shared boundary elements to form a topological complex ^1] for the bounding 
2-manifold without boundary. 

The inherent computational difficulty was to separately create two spline patches, 
each being a manifold with boundary, so that the corresponding boundary curves 
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were identical and could be exactly shared between the patches. In some situations, 
algorithms for fitting spline patches were used successfully. In other cases, patches 
have been slightly enlarged and intersected so as to obtain improved fits. Indeed, 
such intersections are well-defined in pure mathematics, but, again, approximation 
in computation poses subtle variations from that theory, as described in the next 
section on pairwise surface intersection. 

3. Subtleties of pairwise spline surface intersection 

It is well-known that unwanted gaps between spline surfaces or self-intersections 
within intended manifolds often appear as unwanted artifacts of various imple- 
mented intersection algorithms |10| . The mismatch between approximate geometry 
and exact topology has historically caused reliability problems in graphics, CAD, 
and engineering analyses — drawing the attention of both academia and industry. 
The severity of the problem increases with the complexity of the geometric data 
represented, both from high-degree nonlinearity and from the intricate interdepen- 
dence of shape elements that should, but do not, fit together according to the 
specified topological adjacency information. 

The conceptual view of these joining operations is illustrated in Figure|31 with an 
intersection curve^ denoted as c. But this picture of c only exemplifies the idealized, 
exact intersection curve. For practical computations, an approximation of c is often 
created and, in many systems, an intersection curve will be approximated twice; 
once within the parametric domain of one of the intersecting surfaces and then again 
within the parametric domain of the other. These approximations are labeled as 
ci and C2 in Figure |3I The spline functions from [0, 1]^ into IR.^ then also rely upon 
algorithmic evaluation of these approximated intersection curves, as indicated by F 
and G in Figure 13 It is virtually certain that those evaluations will not be exactly 
equal in . 

The mismatch between concept and reality depicted in Figure |3 creates ambigu- 
ity, as the intersection representation is sometimes considered as a unique set, from 
the symbolic topological view, and at other times as two approximating sets, from 
the geometric view. 

4. Specific progress 

To resolve this ambiguity, we are investigating richer representations. We intro- 
duce a neighborhood of the true, but unknown, intersection set. This neighborhood 
is created from newly determined rigorous upper bounds on the error incurred dur- 
ing efficient intersection approximations. To date, it has been convenient to create 
these neighborhoods as tubular neighborhoods |13| . but broader generalizations 
seem to be possible. 

The role of topological equivalence. Before introducing these bounding neigh- 
borhoods, we discuss the meaning of topological equivalence between an object and 
its approximations. We have proposed the use of ambient isotopy for this topolog- 
ical equivalence versus the more traditional equivalence by homeomorphism, as is 
explained more fully in our publications |21I2I'4, 28, 29, 31 . Intuitively, two closed 
curves will not be ambient isotopic if they form different knots. Figure 0] shows 



We focus on the generic case of an intersection curve, although isolated points and co-incident 
areas can also arise, with similar complications. 



COMPUTATIONAL TOPOLOGY FOR REGULAR CLOSED SETS 



5 




Figure 3. Joining operations for geometric objects 



two simple homeomorphic space curves, where the piecewise Hnear (PL) curve is 
an approximation of the smooth curve. However, these curves are not ambient iso- 
topic, because they depict different knots'^, with the smooth curve iUustrating the 
simplest knot, known as the unknot. In the right half of Figure 0] the z coordinates 
of some vertices are specifically indicated to emphasize the knot crossings in 
(All other end points have 2: = 0). All end points of the line segments in the PL 
approximation arc also points on the original curve. Having this knotted curve as 
an approximant to the original unknot would be undesirable as output from a curve 
approximation algorithm, particularly for applications in graphics and engineering 
simulations. Similar pathologies can happen in approximating 2-manifolds, both 
with and without boundary, but results [3 EI El I2H1 12^] summarized here can pre- 
vent these difficulties by appropriately constraining the approximations produced. 

In response to the example of Figure ^ a theorem was published that provided 
for ambient isotopic PL approximations of 1-manifolds jl9|. The proof utilizes "pipe 
surfaces" from classical differential geometry |21J to build an appropriately small 
tubular neighborhood such that if the PL approximant is constrained to lie within 
the constructed neighborhood, then the PL approximant is ambient isotopic to the 
original curve. While the techniques of tubular neighborhoods are well known in 
differential topology, the relevant theorems usually state only the existence of an 
ambient isotopy. Additional work in applied mathematics was needed to comple- 
ment the theorems with constructive formulations from which to obtain effective 
procedures and algorithms. Such extensions were successfully accomplished for 
both non-singular compact, orientable 1-manifolds and 2-manifolds (with or with- 
out boundary) PIH [21123 123 by I-TANGO team members (and their co-authors). 
We used geometric characteristics to compute a specific upper bound on the size of 
a tubular neighborhood and to then specify a particular isotopy. 

^The different knot classifications of Oi and 4™ are indicated in Figure HI 
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Figure 4. Nonequivalent knots 



The importance of this topological equivalence class extends beyond the mani- 
folds described to the geometric models created as compact subsets of 7^(M'^). Con- 
sider that a well-defined topological complex could be created from 2-manifolds 
with boundary if the difficulties along the intersection boundaries could be solved. 
However, the non-uniqueness of the geometric representation of the intersection sets 
seems to pose an intractable problem to creation of a single, well-defined topological 
complex. Our alternative approach is to find a neighborhood within which it can be 
proven that the true intersection curve lies. Since the construction of these bound- 
ing neighborhoods is dependant upon the specific intersection algorithm used, some 
further details of the two intersection algorithms used within the I-TANGO project 
are presented. While the neighborhood construction details will vary with the spe- 
cific intersection algorithm chosen, the following two neighborhood constructions 
were chosen both for their carefully defined error bounds and for their potential for 
generalization. 

Error bounds for topology from Taylor's theorem. First, we present the 
Grandine-Klein (GK) intersection algorithm jl2j. Referring to Figure |2| we note 
that the GK algorithm bases its error bounds on well-established numerical tech- 
niques in differential algebraic equations (DAE). While these DAE techniques pro- 
vide rigorous error bounds, these bounds are expressed within the parameter space 
[0, 1]^, which serves as the domain of the spline functions (indicated as F and G 
in Figure 01 . The code implementing the GK algorithm then has an interface that 
allows the user to specify an upper bound e for this error in parameter space and 
the algorithm provides guarantees for meeting this error bound. However, the typ- 
ical end user is often unaware of the role of this parametric domain, so selection 
of this parametric space error bound has often relied upon heuristics. It would be 
more convenient for the user to be able to specify an error bound within M'^. One 
accomplishment within the I-TANGO project has been to demonstrate a mathe- 
matical relation [221 between the error bounds in M.^ and [0, 1]^, following from a 
straightforward application of Taylor's Theorem in two dimensions [SJ p. 200]. The 
conversion between these error bounds has been implemented in a pre-processing 
interface to the GK algorithm and this new interface has been tested to be reliable, 
efficient and user-friendly. 
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Using the notation from Figure |21 for the sphne function F, Tayfor's Theorem 
provides a bound on the error of F evaluated at a particular point (u, v) versus 
F evaluated at a point (uq^vq), where (w, w) and (uo,vq) are within a sufficiently 
small neighborhood. This sufficiently small neighborhood will have radius given by 
the value in the parametric domain [0, 1]^ which was denoted as e in the previous 
paragraph. Then it follows j22| . with || • || being any convenient vector norm, that 

\\F{u,v)~ F{uo,vo)\\ < eM 

for any AI satisfying 



dF 
ov 



< M, 



for some point on the line segment joining [u,v] and 

For the single spline F, let j{F) be an upper bound for the acceptable error in 
M'^ between the true intersection curve c and one of its approximants F{ci). In 
order to guarantee that this error is sufficiently small, it is sufficient 

eM < 'y{F), 

where an upper bound for Af can be computed by using any standard technique 
for obtaining the maximums of the partials ^ and For G, a similar relation 
between 7(G) and e exists'*. 

Then it is clear that a neighborhood can be defined that contains the true in- 
tersection curve c and both of its approximants. Let F{ci) denote the image of ci 
under F and similarly, let G(c2) denote the image of C2 under G. Let 7V^(p'')(F(ci)) 
be a tubular neighborhood of radius j{F) about F(ci), where ci has been gener- 
ated from the GK intersector to satisfy the inequality presented in the previous 
paragraph. Similarly, define Nj(^q){G{c2))- Then, let 

7V(c)=iV^(;.)(F(ci))UiV^(G)(G(c2)). 

It is clear that N{c) is a neighborhood of c, which contains both of its approx- 
imants, F(ci) and G(c2). However, there is both a theoretical and computational 
limitation to this approach. 

• There is no theoretical guarantee that either approximant is topologically 
equivalent to c, and 

• Any practical computation of N(c) would depend upon an accurate com- 
putation of the set Ny(^p-j{F{ci)) H N^(^g){G{c2)), which is likely to be as 
difficult as the original computation of _F n G. 

While the above bounds are often quite acceptable in practice to compute a rea- 
sonable approximant, further research has been completed into alternate methods 
to give guarantees of topological equivalence within a computationally acceptable 
neighborhood of the intersection set, as reported in the next subsection. 

Integrating error bounds and topology via interval solids. Recent work 
by Sakkalis, Shen and Patrikalakis |3U| emphasized that the numeric input to any 
intersection algorithm has an initial approximation in the co-ordinates used to 
represent points in M.^, leading to their use of interval arithmetic [21]. The basic 
idea behind interval arithmetic is that any operation on a real value v is replaced 



^This error bound assumed that the error due to algorithmic truncation within the numerical 
DAE methods dominated any other computational errors. 
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by an operation of an interval of the form [a, b], where a, 6 G M and a < v < b. The 
resuh of any such interval operation is an interval, which is guaranteed to contain 
the true result of the operation on v. This led naturally to the concept of an interval 
solid and some of its fundamental topological and geometric properties were then 
proven, as summarized below. 

Throughout this section, a box is a rectangular, closed parallelepiped in M.^ with 
positive volume, whose edges are parallel to the co-ordinate axes^. Let be a 
non-empty, compact, connected 2-manifold without boundary. Then the Jordan 
Surface Separation Theorem asserts that the complement of F in M"^ has precisely 
two connected components, Fj, Fq] we may assume that Fj is bounded and Fq is 
unbounded. Let also B = {bj,j G J} be a finite collection of boxes that satisfies 
the following conditions: 

CI: {lnt{bj),j G J} is a cover of F. 

C2: Each member b of B intersects F generically; that is, bOF is a. non-empty 
closed disk that separates b into two (closed) balls, Bj^ and Bj^ , with B^, 
(J5^ ) lying in _F/ U F (Fq U F) , respectively. 
C3: For any bi^bj G B, let 6^ = b., n bj. If Int(fei) n Int(6j) ^ 0, then bij is 
also a box which satisfies C2. 
Notice that condition C2 indicates that every b £ B intersects F in a natural 
way (see Figure EJ. 




Condition C2 Condition C3 



Figure 5. 2D versions of conditions C2 and C3 

The following result summarizes several previously appearing results, where a 
solid is defined to be a non-empty compact, regular closed subset of M"^. 

Theorem 4.1 f |30[ Corollary 2.1, p. 165]). If F is connected and B satisfies Cl- 
C3, then F U Ujgj solid. 

Bisceglio, Peters and Sakkalis [2111231 have recently given sufficient conditions to 
show when the boundary of an interval solid is ambient isotopic to the well-formed 
solid that it is approximating, as described in the following theorem. For a positive 
number (5, define the open set 

F{5) = {xe R^\D{x,F) < S}, 

where D(x^F) — ii\i{d(x^y) \ y G F}, with d being the usual Euclidean metric in 
R3. 



'Enclosures other than boxes are quite possible and this is a subject of active research. 
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Theorem 4.2. Let F be a connected 2-manifold without boundary. For each e > 0, 
there exists S, with < S < p so that whenever a family of boxes B satisfies 
conditions C1-C3, and for each bofB,bisa subset of F(5) (Please see FigureE)) 
then, for S — FUFj and S'^ = SU[Jj^ j bj, the sets F and dS^ are e-isotopic with 
compact support. Hence, they are also ambient isotopic. 

The quoted theorem depends upon results from Bing's book on PL topology |7| 
p. 214], and related literature ^2], as is explained in full [2111221 ■ The proof shows 
that normals to F do not intersect within the constructed tubular neighborhood, 
as is illustrated by the depiction of its planar cross-section in Figure El 




Figure 6. 2D version of proper subset condition 

If the boxes containing the true intersection curve can be made sufficiently small 
so that each such box fits inside F{p), then the resultant intersection neighborhood 
will contain an object that is both close to the true solid and is ambient isotopic 
to it. Considerable success in meeting these constraints has already been achieved 
|24| when two splines intersect transversally, while very recent progress from the I- 
TANGO team for more subtle spline intersection pairs is now under review |23j . As 
a further note on integration, results on root computations from interval arithmetic 
are used to provide estimates of initial starting points for the GK algorithm. 

Work in Progress. The previous discussion emphasizes results from the I- TANGO 
project that have already appeared in the literature, whereas the remarks in this 
section are intended to provide some indication of related results that are expected 
to appear soon. 

Recall that the result of the GK algorithm will be a set of spline patches that 
do not fit together precisely along the approximation of their boundary. Recent 
work by Andersson, Stewart and Zidani uses the Whitney Extension Theorem to 
formulate a conceptual model of how imperfectly fitting patches might be perturbed 
to form a topological complex which is the boundary of a non-empty compact 
element of Tl{R.'^). Furthermore, the proposed process extends the patches under 
Lipschitz mappings so that the resulting element of 7?.(R'^) can be shown to lie 
within rigorous error bounds of the given geometric input data. Our initial results 
on ambient isotopic approximations of 2-manifolds have followed the prevailing 
simplification of considering only those manifolds without boundaries 01211123 j but 
our recent results have been for the more technically challenging cases of 2-manifolds 
with boundaries ^ . Our team is also completing work on spline intersections with 
multiple roots |211, to improve the rather loose error bounds known to-date. 
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Interestingly, the interval solids were not initially intended as a means to in- 
tegrate issues of topological equivalence and approximations in intersection algo- 
rithms. Yet it is a hallmark of our project that these concepts are merging, as 
discussed above. Note, further, that the intervals include both approximation er- 
rors from truncation of numerical processes within intersection algorithms as well 
as the approximations that arise by using a finite set of floating point numbers in 
computation as an approximation of the reals. How floating point approximations 
effect error bounds for intersection computations is another important issue this 
project considers. Attention has focused upon the impact of floating point arith- 
metic on polynomial computations near multiple roots 1161 . This is an important 
issue because the intersection algorithms have to search for roots of the system of 
polynomial equations representing the intersecting surfaces and exact arithmetic 
computations have not been shown to be practical. 

5. Conclusions and future work 

Topology and computer science are finding common interest in the emerging area 
of computational topology. Various branches of pure topology (point-set topology, 
differential topology, low-dimensional topology, . . .) can make important contribu- 
tions to establishing the appropriate theoretical foundations. A fundamentally new 
perspective arises from the role that computational approximations should play 
in the reformulation of central topological concepts. The I-TANGO project is a 
research effort concentrating upon these issues specifically with respect to surface 
intersections within R'^. Considerable progress has already been made, but many 
questions also remain open, as is summarized in this article. 
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